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We measure the top quark mass in dilepton final states of tt events in pp collisions at y/s = 
1.96 TeV, using data corresponding to an integrated luminosity of 9.7 fb _1 at the Fermilab Tevatron 
Collider. The analysis features a comprehensive optimization of the neutrino weighting method to 
minimize the statistical uncertainties. We also improve the calibration of jet energies using the 
calibration determined in tt —» lepton+jets events, which reduces the otherwise limiting systematic 
uncertainty from the jet energy scale. The measured top quark mass is mt = 173.32 ± 1.36(stat) ± 
0.85(syst) GeV. 

PACS numbers: 12.15.Ff, 14.65.Ha 


I. INTRODUCTION 

The discovery of the top quark in 1995 E,a com¬ 
pleted the three quark families of the standard model 
(SM). Since then, the top quark has been one of the focal 
points of the Fermilab Tevatron and of the CERN LHC 
programs. The top quark stands out because of its large 
mass, mt, which is a fundamental parameter in the SM. 
Its Yukawa coupling to the Higgs boson, Y t = \f2m t /v, 
where v is the vacuum expectation value of the Higgs 
field, is close to unity, implying that the top quark may 
play a special role in electroweak symmetry breaking. In 
addition, m t is linked to the W and Higgs boson masses, 
Mw and Mfj, through radiative corrections [sj. Follow¬ 
ing the Higgs boson discovery Si , a precise measure¬ 
ment of nit provides a test of the electroweak sector of 
the SM and information on whether our universe resides 
in a stable or metastable region of that theory Si . The 
short lifetime of the top quark prevents its confinement 
in the strong color field, since top quarks decay before 
lradronizing. This allows a particularly precise study of 
pure quantum chromodynamic (QCD) effects. A com¬ 
parison of the measured m t and the mt extracted from 
cross section measurements Sim may provide a probe 
of higher order and soft QCD corrections to the observed 
mass [ij . 

Assuming the SM branching ratio of t —> Wb « 100%, 
tt decays yield distinct final state categories according 
to the number of charged leptons with high transverse 
momentum (pp) from W boson decays. Dilepton ( 21 , 
I = e or /i) events, such as ee, ep, and p/i, with neu¬ 
trinos from two W —> tv decays, are relatively rare but 
have low background. We present a measurement of mt 


using pp collider data collected with the DO detector at 
the Fermilab Tevatron collider, corresponding to an in¬ 
tegrated luminosity of 9.7 fb _1 , in events with two high- 
Pt electrons or muons of opposite electric charge. Two 
high-pr jets must also be observed, one of which must 
be identified as being consistent with originating from a 
b quark. This analysis is based on our previous dilepton 
measurement (l3| . but with increased integrated lumi¬ 
nosity and multiple optimizations to improve the preci¬ 
sion of mt . We reduce the dominant statistical contribu¬ 
tion to the uncertainty on mt through an optimization 
of the methods for kinematic reconstruction and statis¬ 
tical analysis. Lacking a dijet signature from W -+ qq', 
which is present in tt —> lepton+jets (£+jets) events and 
was used to improve the precision of jet energy calibra¬ 
tion with a W mass constraint [l5j . previous dilepton 
analyses at the Tevatron have reached a sensitivity limit 
imposed by standard jet calibration methods [M E3- 
Pro gres s in calibrating jet energies in the dilepton chan¬ 
nel [t| provides improved cross-checks across different 
channels and a more significant contribution from the 
dilepton channel to the world average mt [H|. For com¬ 
parison, the most recent measurements of mt in the 
dilepton channel from CDF, ATLAS, and CMS are, re- 
spectively, m t = 171.5 ± 1.9(stat) ± 2.5(syst) GeV [19l |. 
m t = 173.79 ± 0.54(stat) ± 1.30(syst) GeV 0, and 
rn t = 172.50 + 0.43(stat) ± 1.46(syst) GeV [2l|. In this 
analysis, we substantially reduce the otherwise dominant 
uncertainty in the jet energy scale by applying the meth¬ 
ods of Ref. [l4j ■ 
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II. DETECTOR AND DATA SAMPLE 

A. Detector 

The DO detector [HI, (HI has a central-tracking sys¬ 
tem, consisting of a silicon microstrip tracker and a cen¬ 
tral fiber tracker, both located within a 1.9 T supercon¬ 
ducting solenoidal magnet, with designs optimized for 
identification of the pp collision vertex and track recon¬ 
struction at pseudorapidities [H[ of \r]\ < 3 and |? 7 | < 2.5, 
respectively. The liquid-argon/uranium calorimeter has 
a central section covering \r]\ < 1 .1, and two end sections 
that extend coverage to \r]\ ~ 4.2, with all three housed 
in separate cryostats. An outer muon system, covering 
| 77 | < 2 , consists of a layer of tracking detectors and scin¬ 
tillation trigger counters in front of 1.8 T iron toroids, 
followed by two similar layers after the toroids. 

B. Object Reconstruction 

We require electrons to satisfy an identification crite¬ 
rion based on boosted decision trees [Hj using calorime¬ 
ter and tracking information. Muons must satisfy re¬ 
quirements that match hits in the muon system to a 
track in the central tracking detector that is required 
to have a small distance of closest approach to the beam 
axis |26j . We require hits in the muon layers inside and 
outside the toroid. Muons and charged hadron momenta 
are measured in the central tracking detector, while elec¬ 
tron, photon ( 7 ), jet, and charged hadron energies are 
measured in the calorimeters. Muons must be isolated 
from jets and from nearby tracks. Electrons and muons 
must have their extrapolated track trajectories isolated 
from calorimeter energy depositions greater than an en¬ 
ergy threshold. Electrons and muons must havepr >15 
GeV, and \p\ < 2.5 and < 2.0, respectively. We recon¬ 
struct jets using an iterative, midpoint-seeded cone algo¬ 
rithm with a cone parameter of lZ CO ne = 0.5 j27|- Jets 
with embedded muons from the decay of 6 -hadrons re¬ 
quire an additional correction to jet energy to account for 
the associated neutrino. A multivariate discriminant [H| 
is used to identify jets that contain a b- hadron (i.e., b 
jets) from a vertex displaced from the interaction point. 
We define the missing transverse momentum (|?r) at¬ 
tributed to the escaping neutrinos as the negative of the 
vector sum of all transverse components of calorimeter 
cell energies, corrected for the measured muon momenta 
and the response of the calorimeter to electrons. We also 
correct for detector response in the jet energy calibra¬ 
tion, as described below. Details of object reconstruction 
are provided in Ref. [29| . 

C. Standard Jet Energy Calibration 

We calibrate the energy of jets to be the energy of 
the particle jets reconstructed using the midpoint algo¬ 


rithm [27]]. We correct for the effects of the calorimeter 
response to particle constituents of jets, energy leaking 
into the cone from particles directed outside it, as well 
as energy deposits outside the cone from particles in¬ 
side it [30(. Charged hadrons have an energy-dependent 
response that is lower than that of electrons and pho¬ 
tons. We therefore apply corrections obtained from 7 +jet 
events to account for the energy dependence of the jet 
response in the central \p\ region. We also apply a rel¬ 
ative 77 -dependent correction obtained from 7 +jet and 
dijet events. We employ the same methods to calibrate 
jet energies independently in the Monte Carlo (MC) sim¬ 
ulation and in data. The MC is used to help study po¬ 
tential biases in the data. We incorporate a correction 
for jets in the MC simulation that accounts for the dif¬ 
ference in single-particle response between data and MC. 
This procedure ensures that the flavor dependence of the 
jet response in data is replicated in MC. In the MC we 
account for multiple pp interactions by correcting the jet 
energy to the particle level of only those particles that 
are directed within the jet cone at particle level. The 
typical systematic uncertainty in the energy calibration 
of each jet in the dilepton sample is 2%. This precision is 
limited by systematic uncertainties of the 7 +jet method 
in the px range of jets in tt events. Details about this 
“standard jet energy scale” calibration can be found in 
Ref. :3(l . We require that jets have px > 20 GeV and 
1 77 1 < 2.5 after calibration, but before applying additional 
corrections from the W -+ qq' constraint in the £+jets 
channel discussed below. 


III. ABSOLUTE JET CALIBRATION FROM A 
W ->■ qq' CONSTRAINT 

As in Ref. [l4|, we apply a multiplicative correction 
factor to the energy of jets in data based on an anal¬ 
ysis of tt— s-f+jets events using the W —> qq' decays 
as a constraint. Application of this factor, 1.0250 ± 
0.0046 (stat) [r|, improves the agreement between MC 
and data and allows us to use its uncertainty to reduce 
the uncertainty on the absolute energy scale by a factor 
of « 4 relative to the standard jet energy scale, while 
retaining its 77 and px dependence. To apply this scale, 
which comes from light-quark jets, to the dilepton sam¬ 
ple, which has b jets, it is important to ensure that the 
variation in the ratio of data over MC jet response be¬ 
tween different flavors be plac ed on an equal footing. The 
standard jet energy scale [301 achieves this on a jet-by-jet 
basis by using single particles in MC jets to correct. This 
ensures that the energies of b jets in dilepton simulated 
samples agree with those of b jets in the dilepton data 
sample at the same level as light-quark jets. Aside from 
fragmentation differences between data and MC which 
are discussed below, this approach justifies the use of the 
£+jets constraint in the dilepton channel. 
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IV. EVENT SELECTION 

The ti candidate events in the ee and pp channels are 
required to pass single-lepton triggers. The full suite 
of triggers is used for selecting ep events. The dilep¬ 
ton event selection before optimization is described in 
Ref. [29|. We optimize the selection based on MC events 
to provide the smallest expected statistical uncertainty in 
m t . We require two isolated leptons with opposite elec¬ 
tric charge. We require at least two jets, where at least 
one of the two jets with highest pr must be identified as 
a b jet. For the ep channel, our selections have an effi¬ 
ciency for tagging b jets of 72%, and a light-quark mistag 
rate of 12% in the central region in p. The same-flavor 
channels employ slightly tighter b tagging requirements 
and thus have a few percent lower efficiency, and 30% 
lower mistag rate. We require events in the pp channel 
to have $t> 40 GeV. This selection is also applied 
to ee events when the dielectron invariant mass is be¬ 
tween 70 and 100 GeV, to reduce the Z —> ee background 
contribution. We define a significance variable, S , 
which measures the likelihood for the observed to be 
a fluctuation from = 0 GeV. We require S > 3.5 (4) 
for the ee ( pp ) channel. We require ep events to have 
Ht > 100 GeV, where Ht is the scalar sum of the pr of 
the two highest-pT jets and of the lepton with highest pr- 
The Ht, b tagging, and ^T-based requirements are opti¬ 
mized to minimize the expected statistical uncertainty on 
mt in each channel. The expected signal-to-background 
(S/B) ratio is « 7 for these channels. These requirements 
yield a 3% improvement in statistical precision in to* rel¬ 
ative to the selections in Ref. [p4|. After implementing 
all these selections, we obtain 340, 115, and 110 events 
in the ep, ee and pp channels, respectively. 


V. MODELING SIGNAL AND BACKGROUND 

The tt events are simulated at 15 mass points over the 
range 130 < mf lG < 200 GeV using the tree level gen¬ 
erator ALPGEN 2.11 [3l| with up to 2 additional light 
partons and PYTHIA 6.409 [32] with modified underlying 
event Tune A for parton showering and hadronization. 
Here, m^ c refers to the input mass in ALPGEN. An ad¬ 
ditional, larger sample is generated at mf lG = 172.5 GeV 
to study systematic uncertainties. We normalize the tt 
production cross section to a t t = 7.24 ± 0.23 pb [33|, 
which is calculated at next-to-next-to-leading order with 
a next-to-next-to-leading logarithm soft gluon resumma¬ 
tion. The main backgrounds arise from three sources: 
Z/ 7 * —> i + i ~, diboson (WW, WZ, and ZZ) processes, 
and instrumental effects. We model the Z /y* background 
using ALPGEN with up to 2 light partons and PYTHIA for 
showering and hadronization. We employ PYTHIA for 
the diboson background. The instrumental background 
arises from VF+jets, multijet, or ^-f-jets tt events where 
one or two jets are either mis-identified as electrons, or 
they contain a hadron decaying to a non-isolated lepton 


that passes our selection. This background is estimated 
from data as in Ref. [2§|. We apply a full detector simu¬ 
lation based on GEANT 3.14 [3J] for all simulated events. 
The objects reconstructed in simulation are smeared to 
ensure that their resolutions reflect those in data. Scale 
factors in object efficiencies are applied to improve agree¬ 
ment between data and MC. 


VI. KINEMATIC RECONSTRUCTION 


A. Neutrino Weighting 


The presence of two neutrinos in the ti decay makes it 
impossible to fully constrain the kinematics and thus ex¬ 
tract a unique mt measurement from each event. Given 
the measured momenta of leptons, jets and $t, the avail¬ 
able constraints from Mw, and the condition mt = 
we are missing one constraint to provide full tt recon¬ 
struction in dilepton events. We integrate over the phase 
space of neutrino rapidities for chosen values of hypoth¬ 
esized m t (m 1 ^) HH, and compare ^|? lc , the vector sum 
of neutrino momentum solutions at each chosen point 
of phase space, to the observed to determine a 
“weight” lo characterizing the level of agreement: 


N 


inn 


exp -- 


i=l j=x,y 


/77'calc z7 , obs\2 
VFT 7 - fi $Tj ) 


(1) 


where i runs over all neutrino solutions for any two pos¬ 
sible jet-lepton assignments in the ti final state (up to 
N = 8 ), j stands for the two orthogonal coordinates in 
the transverse plane (x and y ), and is the resolu¬ 
tion of the x or y component of unclustered The 
unclustered $t ($t) is the magnitude of the vector sum 
of all energy depositions in the calorimeter that were 
not included in lepton or jet reconstruction. The quan¬ 
tity is parameterized to describe the difference be¬ 
tween the neutrino momentum solutions and the mea¬ 
sured in the weight calculation. We perform this 
calculation over a range of rn\. integrating u; over the 
neutrino phase space, to yield a distribution of ui(mt) 
versus my. Prior studies [361 ] have shown that the first 
two moments (p u ,cr u ) of this distribution extract most 
of the information about mt- The analysis of Ref. fl4[ 
used the range of mj 1 values between 80 and 330 GeV in 
1 GeV step and a resolution of 7 GeV in the weight 
calculation. The new optimized determination of these 
parameters is briefly summarized below. 


B. Optimization of Weight Calculation Parameters 

After applying the methods described above to im¬ 
prove the jet energy calibration, the statistical contribu¬ 
tion is the dominant source of measurement uncertainty 
on mt hr the dilepton channel. We therefore examine the 
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parameters used for the kinematic reconstruction of tt 
events and for the maximum likelihood fit to reduce the 
expected statistical uncertainty. At each step, we verify 
through MC simulations that the optimization does not 
increase the systematic uncertainty. 

All neutrino solutions and jet assignments yield mass 
estimators such as p^ that are correlated with rrit . How¬ 
ever, the correlation is substantially greater, and p u val¬ 
ues are less biased, when the correct jet assignments and 
solutions of neutrino momenta are chosen. Since now 
rrit has been measured with high precision (l8j , we can 
optimize the range of mj 1 based on known values of m*. 
Considering a wide range in m causes incorrect config¬ 
urations to overwhelm the correct configuration, thereby 
worsening the mass resolution. Likewise, scanning over 
too narrow a range biases the background and worsens 
the mass sensitivity by causing tt and background distri¬ 
butions to be similar. Examination of a two-dimensional 
grid of upper and lower limits of the mass range yields the 
optimal range of mj 1 = 115 to 220 GeV in 1 GeV steps. 
The value of also has a noticeable impact on the 
expected precision of the analysis. Properly accounting 
for the full difference between the calculated and mea¬ 
sured j including effects such as mismatches between 
the true neutrino pt and that taken from finite binning 
of the neutrino rapidity distributions, or the difference 
in the number of jets used to derive the calculated and 
measured requires that we optimize the value of the 
resolution parameter. We scan for a wide range 
from 7 to 100 GeV and find 25 GeV is optimal. Com¬ 
bined, these optimizations improve the expected com¬ 
bined statistical uncertainty on rrit by 11% compared to 
the parameters in Ref (l3 |. 


C. Efficiency of Kinematic Reconstruction and 
Event Yields 


Events used in the analysis must have at least one pair 
of neutrino solutions for at least one mj 1 value. The 
efficiency for this kinematic reconstruction is over 99% 
for tt events, and 91% to 98% for the background. In 
the final sample, a total of 336, 113, and 109 events in 
the ep, ee, and pp channels, respectively, pass the kine¬ 
matic reconstruction. The expected sum of tt and back¬ 
ground yields and their corresponding asymmetric total 
uncertainties (stat©syst) are 298.106.5i}^6> an d 
103.5jl9'i events for the ep, ee, and pp channels, respec¬ 
tively. The distributions of the mass estimator p u in a 
preselected sample, omitting requirements on b tagging, 
|?t, significance, and Ht, are shown in Fig. Uta). The 
tt component is evident in the preselected data. The mass 
dependence of the p u distribution is given in Fig[ljb) for 
three mass points with all selections applied. 



120 140 160 180 200 220 

9 w [GeV] 



FIG. 1. The distribution in the mass estimator, p w , for the 
combination of the ee, ep, and pp channels for (a) the prese¬ 
lected sample and (b) the final event sample. The MC events 
are normalized separately to the number of observed events 
in data in each channel. The ratios show the total number of 
observed events divided by the number of expected events in 
a given bin of p w for mp' ic = 172.5 GeV. The band of system¬ 
atic uncertainty is shown as the shaded area in the ratio plots, 
which includes contributions from the dominant sources: jet 
energy scale, lepton identification, lepton momentum scale, 
luminosity, b quark modeling, initial and final state radiation, 
color reconnection, as well as hadronization and higher-order 
QCD effects for tt events. 


VII. EXTRACTING THE TOP QUARK MASS 
A. Maximum Likelihood 

We perform a binned maximum likelihood fit to the 
extracted moment distributions [p u ,<7 u \ in data. Ex¬ 
pected probability densities are calculated using the MC 
samples for each of the 16 rrit points, yielding a two- 
dimensional probability density hs(pu, o©|to^ c ) distri¬ 
bution parametrized by rn t . Background samples are 
used to construct a background template for each chan¬ 
nel, hBiPun&uj), with each background contributing ac¬ 
cording to its expected yield. Bins in signal templates 
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TABLE I. Slopes, offsets, and pull widths of the mt calibra¬ 
tion and the expected statistical uncertainties in the mass 
(am t ) for the ee, efi , and /r/r channels, and their combination. 



Slope 

Offset [GeVJ 

Pull width 

<J mt [GeVJ 

ee 

0.984 ±0.004 

0.671 ±0.043 

0.994 

2.98 

e/i 

0.986 ±0.006 

0.548 ±0.065 

0.998 

1.72 


0.989 ±0.010 

0.717 ±0.103 

1.004 

3.31 

2t 

0.988 ±0.006 

0.617 ±0.063 

0.995 

1.35 


with no events are given a weighted value correspond¬ 
ing to a single signal MC event to ensure that the log of 
likelihood is not infinite. The likelihood is given by: 

N | tls, tlB , Ult ) = 

n N n s ■ hs{li ui ,<Tui I m t ) + n B • 

n S + TIB ’ 

l=L 

where N is the number of observed events in data, ns is 
the expected number of tt events (for m t = 172.5 GeV), 
and ub is the expected total number of background 
events. We fit (— In C) versus rn^ c to a parabola in a 
window of inf LC that is iteratively varied until a stable 
minimum is found. We take the minimum of the final 
parabola to be the fitted top quark mass, m® 1 . The un¬ 
certainty on the fitted mass is obtained by considering 
the m^ c range over which the fit function increases by 
0.5 units in (— ln£) above this minimum. Using pseudo¬ 
experiments, we optimize the template binning of each 
channel separately in a two-dimensional grid that lets 
and cr u bin sizes vary independently. Finer binning in 
and (Tu, especially for the efi channel, improves the ex¬ 
pected statistical precision in m® 1 by 5%. The fitted mass 
window is optimized to ±15 GeV for all channels. Taking 
all the optimizations together, including event selection, 
weight calculation, and maximum likelihood fitting, the 
statistical sensitivity of this analysis is improved relative 
to Ref. 0 by 20% beyond the 35% gain expected from 
increased integrated luminosity. 

B. Ensemble Testing and Data Results 

We obtain a linear relationship between to® 1 and mf LC 
by performing randomized pseudo-experiments using all 
signal mass points. The numbers of signal and back¬ 
ground events in the pseudo-experiments are allowed to 
fluctuate within their Poisson uncertainties around their 
expected values. We require that the total number of 
events matches that observed in data. To minimize the 
effect of statistical fluctuations on our systematic uncer¬ 
tainties, we optimize the number of pseudo-experiments 
by dividing the MC sample into five subsamples, and 
measure systematic uncertainties with each subsample. 
We calculate the RMS of the five uncertainties, aver¬ 
age over all systematic effects, and divide by \/5 to esti¬ 
mate the statistical component of systematic uncertain¬ 


ties. The average RMS decreases until we oversample, or 
reuse, the tt MC events by roughly a factor of three. This 
corresponds to 3000 pseudo-experiments. We perform a 
linear fit of to® 1 versus to^ c to obtain a calibration slope 
and offset for to® 1 using 3000 pseudo-experiments: 

raf 1 = Slope • (m f MC — 170) ± Offset ± 170. (3) 

We account for oversampling by increasing the statisti¬ 
cal uncertainties at each mass point by the appropriate 
oversampling factor. Likewise, we compute the pull, or 
the ratio of to® 1 —to^ 10 over the average estimated un¬ 
certainty at each mass point. The slopes of to® 1 ver¬ 
sus to^ c are close to 1, and pull widths are consistent 
with unity, as shown in Tabic |T| We calculate the final 
to* by correcting to® 1 from a given measurement by the 
slope and offset. We correct the statistical uncertainty 
using the slope and the pull width. The expected cor¬ 
rected statistical uncertainties for each channel are given 
in Table Q] In data, we obtain corrected, fitted mt val¬ 
ues of to 4 = 171.86 ± 1.71(stat), 173.99 ± 3.04(stat), and 
178.58 ± 3.56(stat) GeV for the e/it, ee, and /i/i channels 
respectively, and m t = 173.32 ± 1.36(stat) GeV for the 
combined channels. 


VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties summarized in Table ITT1 arise 
from jet energy calibration, object reconstruction, model¬ 
ing of tt and background events, and the mass-extraction 
method. The energies of jets are shifted up and down by 
the uncertainty on the absolute energy scale, which is 
taken from £±jets events, thereby providing shifts in mt- 
This scale is appropriate for light-quark jets, which, af¬ 
ter correcting for jet flavors to improve the agreement 
between data and MC, have different kinematic distribu¬ 
tions than b jets from tt decays. We calculate a resid¬ 
ual uncertainty due to the kinematic differences between 
the Cfjets calibration sample and dilepton sample of b 
jets. We use separate up and down estimates to extract 
the energy- and ^-dependent shifts in m* based on un¬ 
certainties in the standard jet energy scale relative to 
their average value in the f±jets calibration sample. We 
cross-check this with an alternative method that applies 
shifted light-quark jet energy scales to b jets in the £±jets 
channel [l5| . These methods agree, and thereby validate 
the use of the £±jets scale as a jet calibration. We also 
cross-check using a jet-energy-dependent linear parame¬ 
terization of the residual jet energy scale as in Ref. Il5j . 
obtaining results that do not exceed our estimate of un¬ 
certainties from the jet energy scale. To estimate the 
uncertainty corresponding to possible differences in the 
flavor dependence of the MC scale relative to data, we 
change the single-particle responses up and down by their 
uncertainties and obtain the shift in mt- To estimate 
the possible dependence on the b quark fragmentation in 
the MC, we replace the PYTHIA b quark fragmentation 












TABLE II. Systematic uncertainties on mt for the combined 
dilepton measurement using 9.7 fb -1 of integrated luminosity. 
For symmetrized uncertainties, the “±” symbol indicates that 
the corresponding systematic parameters in MC are positively 
correlated with mt in data, and the “+” symbol indicates an 
anticorrelation. The uncertainties shown as + or — only are 
computed by comparing a standard choice with an alternate, 
but are symmetrized in calculating the total uncertainty. 


Source 

<J mt [GeV] 

Jet energy calibration 


Absolute scale 

+0.47 

Flavor dependence 

TO.27 

Residual scale 

+0.36 

-0.35 

b quark fragmentation 

+0.10 

Object reconstruction 


Trigger 

-0.06 

Electron px resolution 

±0.01 

Muon pt resolution 

+0.03 

Electron energy scale 

±0.01 

Muon p T scale 

±0.01 

Jet resolution 

+0.12 

Jet identification 

+0.03 

b tagging 

+0.19 

Signal modeling 


Higher-order effects 

-0.33 

ISR/FSR 

±0.15 

Pr(tt) 

-0.07 

Hadronization 

-0.11 

Color reconnection 

-0.22 

Multiple pp interactions 

-0.06 

PDF uncertainty 

±0.08 

Background modeling 


Signal fraction 

±0.01 

Heavy-flavor scale factor 

±0.04 

Method 


Template statistics 

±0.18 

Calibration 

±0.07 

Total systematic uncertainty 

±0.85 


function with the Bowler scheme [37}, and compare mt 
with the Bowler free parameters tuned to LEP (ALEPH, 
OPAL, and DELPHI) or SLD data H|. 

The systematic uncertainty due to the trigger efficiency 
is estimated by applying the ratio of single lepton trigger 
efficiency parameterization in data divided by the MC pa¬ 
rameterization to the ee and pp channels. The uncertain¬ 
ties in the modeling of the energy and momentum reso¬ 
lutions of electrons, muons, and jets are applied indepen¬ 
dently of each other, and the shifts in to* are extracted as 
uncertainties on mt- Lepton energy or momentum scales 
and their uncertainties are extracted from Z —> 2£ events 
in data. An additional uncertainty is estimated for jet 
identification by shifting the jet identification efficiency 
within its uncertainty in MC samples to estimate their 
effect on m t . The uncertainty from modeling b tagging 
is evaluated by changing within their uncertainties the 


corrections that account for the agreement between data 
and MC in b tagging efficiency. 

Higher-order virtual corrections to m t are absent in the 
ALPGEN used to generate our standard tt samples. We 
therefore compare an ensemble of pseudo-experiments us¬ 
ing MC@NLO 3.4 [H] tt events with one using ALPGEN 
events, where both employ HERWIG 6.510 [40] for mod¬ 
eling of hadronization. To evaluate the uncertainty as¬ 
sociated with the modeling of initial and final-state ra¬ 
diation (ISR/FSR), we compare ALPGEN+PYTHIA with 
the renormalization and factorization scale changed up 
and down by a factor of 1.5 [lj|. The ^+jets analy¬ 
sis exhibits a discrepancy in the shape of the pt dis¬ 
tribution of the tt system, which, although the dilep¬ 
ton statistics are limited, may be present in the dilepton 
sample. We evaluate the uncertainty in the modeling 
of the tt pt distribution by reweighting MC events to 
make them match the data. The observed shift in mt 
is taken as the uncertainty. Since the hadronization in 
our standard tt sample is modeled with PYTHIA, we es¬ 
timate a hadronization uncertainty on m* by performing 
pseudo-experiments using an ALPGEN+HERWIG sample. 
We evaluate the effect of color reconnection by comparing 
mt measurements in ALPGEN+PYTHIA samples with two 
PYTHIA tunes: the Perugia2011 tune that incorporates an 
explicit color-reconnection scheme, and the Perugia2011 
NOCR tune that does not [4l|. Data and MC may have 
different distributions in instantaneous luminosity after 
event selection. This uncertainty due to multiple pp in¬ 
teractions is estimated by reweighting the distribution 
of instantaneous luminosity to make MC agree with the 
data for respective data-taking epochs, and then take the 
shift in mt with respect to the default value. The uncer¬ 
tainty due to the proton structure is obtained from the 
20 sets of CTEQ6L1 parton distribution functions (PDF) 
reweighted to CTEQ6M, where the deviations in nu for 
the 20 eigenvectors sets are added in quadrature |42 |. 

We estimate the effect of the uncertainty on the frac¬ 
tion of signal or background by changing the expected tt 
event yields (ns) up and down and the expected back¬ 
ground yields (ng) down and up within their total uncer¬ 
tainties. The heavy-flavor scale factor, which is applied 
to the Z —> 2£ cross section to correct the heavy-flavor 
content, is also changed up and down within its uncer¬ 
tainty to estimate its systematic effect on mt- 

Our templates are constructed from MC samples for 
tt, Z —I 2£, and diboson backgrounds, as well as data 
samples for instrumental background, yielding statistical 
uncertainties on their bin contents. We use Poisson dis¬ 
tributions to modify bin contents within their statistical 
uncertainties to obtain 1000 new templates. We measure 
mt in data using these templates, and the RMS of the 
measured top quark mass is taken as its uncertainty. Our 
method of mt extraction relies on the correction of the 
fitted mt to the input MC mass. The uncertainties from 
this calibration are applied to provide the uncertainty 
in mt- The uncertainty is reduced substantially from 
Ref. [hi ] due primarily to the reduction in the uncer- 
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tainty in jet energy calibration and the optimizations for 
improvements in statistical uncertainty. Larger MC sam¬ 
ples also contribute by lowering statistical fluctuations on 
systematic uncertainties, or reducing statistically limited 
systematic uncertainties. 

IX. CONCLUSIONS 

We have measured the top quark mass in the combined 
dilepton channels (e/Li, ee, /r/i): 

m t = 173.32 ± 1.36(stat) ± 0.85(syst) GeV 
= 173.32 ± 1.60 GeV. 

This measurement is consistent with the current aver¬ 
age value of mt HI! . Our measurement is the most pre¬ 
cise dilepton result from the Tevatron, and is competi¬ 
tive with the most recent LHC dilepton measurements. 
The systematic uncertainty of 0.49% is the smallest of all 
dilepton measurements. 
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